Charged colloids at low ionic strength: macro- or microphase separation? 
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Phase separation in charged systems may involve the replacement of critical points by microphase 
separated states, or charge-density-wave states. A density functional theory for highly charged 
colloids at low ionic strength is developed to examine this possibility. It is found that the lower 
critical solution point is most susceptible to microphase separation. Moreover the tendency can be 
quantified, and related to the importance of small ion entropy in suppressing phase separation at 
low added salt. The theory also gives insights into the colloid structure factor in these systems. 
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There has been much interest recently in statistical 
physics in charged soft matter systems. Whilst much of 
this is biologically inspired (eg DNA condensation |E|), 
there has also been a long standing controversy in the 
colloid science community over the anomalous behaviour 
of charge stabilised colloidal suspensions at low ionic 
strengths 0. Recently Hf§, it has been suggested 
that the anomalies in these systems may be understood 
in terms of a miscibility gap which is the analogue of 
the vapour-liquid coexistence in the restricted primitive 
model (RPM) of 1:1 electrolytes ||. Arguably a theo- 
retical consensus is emerging, although there remain a 
number of competing theories Q . 

In all these examples, the crucial role of the counterions 
should not be underestimated. For the case of charged 
colloidal suspensions, for example, the theories show that 
the overall phase stability is almost entirely due to the 
entropy of the counterions |3|4|]. The same effect can 
be said to underpin the solubility of many water-soluble 
polymers 0. The basic point is that bulk phase separa- 
tion in a charged system must be into electrically neutral 
phases. If this involves significant fractionation of coun- 
terions, an entropic penalty will be incurred which tends 
to suppress phase separation. 

Clearly though, if bulk phase separation is suppressed, 
a possibility still exists to undergo microphase separa- 
tion, where electroneutrality can be broken locally. Crit- 
ical points are particularly susceptible to this, as first 
shown by Nabutovskii, Nemov and Peisakhovich (NNP) 
using a Landau- Ginzburg theory ||. Consider density 
fluctuations at wavevector q. At q — * 0, fluctuations 
are restricted to elecrically neutral combinations, but for 
q > fluctuations can violate electroneutrality increas- 
ingly easily. Thus one might expect some softening of the 
modes. Indeed, if the only terms to 0(q 2 ) come from the 
long range Coulomb law, the analysis below implies that 
all partial structure factors have a minimum at q = 0. 
Since the q = partial structure factors diverge as one 
approaches a critical point, this suggests there must ex- 
ist regions around critical points where a divergence at 
q > occurs first, indicative that the critical behaviour 
is preempted by microphase separation. However, there 



are often other terms arising at 0(q 2 ) from elsewhere 
which destroy the phenomenon. A closely analogous mi- 
crophase separation for polyelectrolytes in poor solvents 
has also been examined but in the present study mi- 
crophase separation is driven purely by electrostatic ef- 
fects. 

Let us start by constructing a simplified but physically 
motivated model for the anomalous behaviour in charged 
colloidal suspensions. Consider the macroions as spheres 
of charge Z , diameter 2a, and number density pm (vol- 
ume fraction <j) = Ana 3 pm/3). There are small counte- 
rions and coions at number densities p_ and p+ respec- 
tively. The solvent is a dielectric continuum. Without 
loss of generality, I suppose the small ions are univalent, 
and there is only one species of counterion Since 
the coions come from added salt, it will be convenient to 
write p+ = ps. Overall, the system is electrically neu- 
tral and p- — ps + ZpM, but Sp± will be retained for 
fluctuations. 

Each macroion polarises the surrounding electrolyte 
solution, and becomes surrounded by a 'double layer'. 
It has been shown by many workers P, fil]|l2| that the 
self energy of the macroion with its double layer, in 
Debye-Hiickel theory, is {Z 2 l B kT/2a) x h(na). In this 
Ib = e 2 /ekT is the Bjerrum length, the function h(x) = 
1/(1 + x) and the Debye screening length, is 
given by k 2 — 8itIbPi where 2pi = p_+p + = ZpM+2ps is 
(twice) the ionic strength. This energy has a well known 
interpretation: it corresponds exactly to a spherical ca- 
pacitor, charged ±Ze, with one plate at the macroion 
surface and the second a distance kT 1 away [H . 

The simplest model free energy based on this is 

F/VkT = p s logps + (ps + z Pm) log(ps + Zp M ) 

+PMlogpM + Pm (Z 2 l B /2a) h(na). (1) 

The first three terms are the ideal terms, and the last is 
the self energy of macroions at number density pm ■ The 
most important omission from this is the contribution 
from the macroion-macroion interactions. Whilst this 
plays a significant role in structuring the macroions, it 
has been demonstrated elsewhere Q that it is less signif- 
icant than the self energy, as regards the appearance of 
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a miscibility gap. Moreover, by leaving this contribution 
out of the theory, we will see quite clearly how structure 
can develop in the system in the absence of pair interac- 
tions. 

A typical phase diagram corresponding to the above 
free energy is shown in Fig. (jl for Z = 10 3 and 2a — 
100 nm. It comprises a simple miscibility gap, limited 
above and below by critical solution points as the salt 
concentration is varied. The gap occurs at very low 
ionic strengths, and only appears if the charge on the 
macroions is sufficiently high (ZI-q/cl ^ 13.4 for 2o in the 
range 10-1000 nm). 

Now, the NNP scenario could occur at either crit- 
ical point. To examine this therefore, I construct 
a density functional theory to correspond to the free 
energy introduced above. The ideal terms become 
/cT/d 3 rpi(r)logpj(r) (i = +,-,M), and I intro- 
duce the ansatz that the self energy generalises in 
the obvious way to Jd 3 r pm(i") /|/ lf (r) where /^ lf = 
(Z 2 l-QkT '/2a)h(na) is the self energy per particle evalu- 
ated using the local ionic strength at the particle centre, 
k 2 = 8itIbPi(y)- Finally, an electrostatic contrbution has 
to be added: l B kT Jd 3 rd 3 r'pz(r)pz{r')/\r — r'|, where 
pz(r) = Zpwi(r) + p+(r) — p~(r) is the local charge den- 
sity. 

To examine the stability of the system against mi- 
crophase separation, expand the above density functional 
about the homogeneous state to quadratic order. For 
fluctuations at a wavevector q this results in 
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The functions are h\{x) = (1 + 3x)/(x 3 (l + x) 3 ) and 
h2{x) = l/(x(l + x) 2 ). From this the macroion structure 
factor S(q) = (\3pM{q)\ 2 ) is extracted in the standard 
way The behaviour of S(q) is examined as a func- 
tion of 4> and ps, looking for the unstable regions in the 
(0, ps)-plane where l/S(q) < 0. 

At q = the spinodal instability region corresponding 
to the free energy in Eq. (|l|) is recovered. For q > 
the region of instability always expands. This is because 
the (/-dependence arises solely from the long range elec- 
trostatic term, thus, as alluded to above, S(q) always 
has a minimum at q = 0. But a clear difficulty emerges 
when the behaviour for large q is examined, since the 
instability region expands to fill the entire plane; there 
is no effective penalty against microphase separation at 
vanishingly small wavelengths. 
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FIG. 1. Miscibility gap at low ionic strength for model free 
energy, Eq. ([j} in text, at Z = 10 3 and 2a = 100 nm. Also 
marked are the lines where the counterions and added salt 
contribute equally to the ionic strength (ZpM = 2ps), and 
where the macroion self energy function varies most rapidly 
(x = Ka = 1). 

In fact the ansatz for the self energy term has omit- 
ted an obvious but crucial effect, namely one would not 
expect a macroion to be sensitive to variations in the 
local ionic strength over distances much smaller than its 
size. Motivated by weighted local density theories for liq- 
uids Jl6| , I therefore introduce an additional smoothing 
ansatz. It turns out that the precise form does not mat- 
ter greatly, for instance one can replace the local ionic 
strength by pj(r) = Jd 3 r'w(\r — r'\)pi(r') where w(r) 
is a smoothing kernel of range a | p7[ , but equally one 
might smooth k or /j^ lf . All forms result in the appear- 
ance of extra multiplicative factors, wi(qa) and W2(qa), 
in the last two terms of Eq. The Wi(qa) are related 
to the Fourier transform of w(r), and satisfy Wi — > 1 
at q — > 0, Wi — > at q — > oo. (More generally I ex- 
pect hi(na,qa) such that hi — > for qa — > oo.) With 
this ansatz, progress can be made without developing a 
detailed theory by investigating various possibilities for 
Wi. The results reported below have been carried out 
assuming w± = W2 = cxp(—aq 2 a 2 ) with a a numerical 
prefactor of order unity. Very similar results are obtained 
for Wi = 1/(1 + aq 2 a 2 ). 

Typical results from this modified theory are shown in 
Figs. U and |], for a = 1. In Fig. ||, the spinodal insta- 
bility at q = is recovered as before. For q > 0, the 
instability region is again expanded in the vicinity of the 
lower critical point, but is now reduced in the vicinity of 
the upper critical point. For qa ^> 1 the instability dis- 
appears completely, since the self energy which drives the 
instability is now insensitive to short wavelength fluctu- 
ations. Below the lower critical point, therefore, there is 
a region (delimited by the heavy dashed line in Fig. |^) 
where S(q) diverges at some q* > 0, corresponding to the 
appearance of microphases. 
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FIG. 2. Phase diagram in double logarithmic plot, showing 
as dashed lines the contours of divergence of S(q), from q — 
(identical to original spinodal), through ga = 0.05, 0.1, 0.2 
(light dashed lines) to q = q* (heavy dashed line) where the 
maximum in S(q) diverges (typically q*a < 0.5). The struc- 
ture factors corresponding to the marked points are given in 
Fig. |. 

These results are reflected in the macroion structure 
factors, two examples of which are shown in Fig. |^. Near 
the upper critical point the structure factor turns up to 
a maximum as q — > 0, developing a divergence at q = 
as the (mean field) critical point is approached. Near the 
lower critical point a peak appears in S(q) at q* > 0. The 
peak diverges as one approaches the boundary of the mi- 
crophase instability region. Note that the appearance of 
a peak in the macroion structure factor S(q) is unusual 
because there are no direct macroion interactions in the 
theory as constituted above. This shows how structure 
can arise in a charged system independent of the exis- 
tence of (effective) pair interactions. 

To relate the results to possible experiments, one 
should of course investigate the significant contribution 
to the structure factor from the omitted macroion inter- 
actions. Elsewhere it is argued that, at these low ionic 
strengths, the macroions are effectively a one-component 
plasma (OCP) in the strong coupling limit Q, whose 
structure factor resembles that of hard sphere (HS) fluid 
close to the freezing transition |Q. A rescaled S(q) 
for HS at freezing is compared with the present calcu- 
lations in the inset in Fig. ||. The additional structure 
arising from the self energy theory above appears to lie 
well within the first peak, and moreover the amplitudes 
match quite closely. Thus the appearance of a maximum 
at q = or a peak at very low q may well be experi- 
mentally observable. This may be the explanation of the 
anomalously large S(q — > 0) reported recently for col- 
loidal suspensions at low ionic strengths Jl9| . 

Let us turn to the effect of a. Recall that a is a 
measure of the degree of smoothing: the range of the 
smoothing kernel is ~ ay/a. All the results discussed 
above were at a = 1. If a < 0.338, microphases ap- 



pear at the upper critical point too. On the other hand, 
if a > 1.961 microphase separation at the lower critical 
point disappears. One would expect increasing a to sup- 
press microphase separation, since greater smoothing is 
bound to reduce q* , but the difference between the two 
critical points is suggestive. These critical values of a 
are more general than the assumed form of Wi(q), since 
they only depend on the q 2 coefficient in the expansion 
Wi(q) = 1 — aq 2 a 2 + 0(q 4 ) p0| . They are a measure of 
the susceptibility of the critical point to replacement by 
microphases (if a could be treated as a control variable, 
the critical values would correspond to Lifshitz points in 
the phase diagram). 
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FIG. 3. Structure factor in vicinity of miscibility gap. The 
state points are those marked in Fig. 0. The dashed line is the 
structure factor expected from the interparticle repulsions at 
a scaled volume fraction around freezing (all normalised such 
that S(q) — - > 1 at q — > oo). 

This calculation sheds light on the reason for a closed 
loop miscibility gap. As ps is increased, an upper crit- 
ical point is expected since the self energy ceases to be 
strongly state point dependent for na 3> 1 or ps ZpM, 
where it simply shifts the macroion chemical potential 
(see lines in Fig. Q). As ps — > though, it is perhaps 
unexpected to encounter a second critical point. Its ap- 
pearance appears to be connected to the small ion en- 
tropy effect discussed in the introduction. The evidence 
for this is twofold. Firstly, as already commented upon, 
the lower critical point is more susceptible to microphase 
separation. This is in accord with the idea that phase 
separation is suppressed by small ion entropy, since in 
microphase separation, the system gains entropy by dis- 
tributing the small ions more uniformly than would be 
allowed if strict electroneutrality had to be satisfied at 
each point. 

The second piece of evidence concerns the rate at which 
the Donnan potential difference Aip vanishes as one ap- 
proaches the critical point. Recall that Alp arises be- 



3 



cause the interface can acquire a dipole moment den- 
sity (per unit area). In the present system, a dipole 
moment density appears to arise because the jump in 
small ion densities is spread out more broadly than the 
jump in macroion densities (although this remains to be 
confirmed with a detailed calculation Remarkably, 
one can calculate Atp without detailed knowledge of the 
structure of the interface Q. I find that Aip vanishes 
as A</>/(/> cr it as the critical points are approached, with a 
constant of proportionality « 8.00 for the upper critical 
solution point, and 16.9 for the lower one. This again 
indicates the growing importance of small ion entropy 
(broadening the jump in small ion densities) as the lower 
critical point is approached. 

Note that Aip is an order parameter which strictly 
vanishes in symmetric models such as the RPM. Apart 
from general remarks by Nabutovskii and Nemov |2^] , the 
critical behaviour of asymmetric primitive models seems 
to have received much less attention than the RPM ||, 
and there may be interesting effects connected to a non- 
vanishing Atp. 

I thank C. Holmes and M. E. Cates for useful discus- 
sions, and P. Schurtenberger for sending data prior to 
publication. 
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